## Construction of F16 in two steps

from PyM import *

F2 = Zn(2)

[F4,x] = extension(F2,[1,1,1],'x')

# The polynomial Y**2 + x*Y + 1 is irred over F4 (has no roots in this field)
[F16,y] = extension(F4,[1,x,1],'y')


show(order(y))

show(order(y+1))

show(geometric_series(y+1,cardinal(F16)))